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Dynamics of dissipation of a local phonon distribution to the substrate is a key issue in friction 
between sliding surfaces as well as in boundary lubrication. We consider a model system consisting 
of an excited nano-particle which is weakly coupled with a substrate. Using three different methods 
we solve the dynamics of energy dissipation for different types of coupling between the nano-particle 
and the substrate, where different types of dimensionality and phonon densities of states were also 
considered for the substrate. In this paper, we present our analysis of transient properties of energy 
dissipation via phonon discharge in the microscopic level towards the substrate. Our theoretical 
analysis can be extended to treat realistic lubricant molecules or asperities, and also substrates 
with more complex densities of states. We found that the decay rate of the nano-particle phonons 
increases as the square of the interaction constant in the harmonic approximation. 



I. INTRODUCTION 

Friction between two surfaces in relative mo- 
tion involves many interesting and complex phenom- 
ena induced by the long- and short-range forces, 
such as adhesion, wetting, atom exchange, bond 
breaking and bond formation, elastic and plastic 
deformatior^^iS&i ^ 2 .^.. In general, non- 
equilibrium phonons are generated in the expense of 
damped mechanical energ y 15 i 16 ' 17 i 18 i 19 i 20 . Dissipation of 
this excess energy is one of the important issues in dry- 
sliding friction and lubricatio n 21 i 22 i 23 i 24 . Normally, the 
dissipation of mechanical energy is resulted in heating 
of parts in relative motion. Sometimes, it gives rise to 
wear and failure due to overheating. In general, signifi- 
cant amounts of resources (energy and material) are lost 
in the course of friction. One of the prime goals of tri- 
bology is to minimize energy dissipation through lubrica- 
tion. Recently, several works have attempted to develop 
surfaces with superlow friction coefficients. 

The objective of this work is to develop an under- 
standing of phononic energy dissipation during sliding 
friction, especially to shed some light on the dynamics 
of discharge of excited phonons of a nano-particle (rep- 
resenting a lubricant molecule or an asperity) deep into 
the substrate. This problem has many aspects and the 
solution will depend on a variety of physical parameters 
which can be grouped into major categories, such as in- 
ternal degrees of freedom of the nano-particle, density of 
substrate phonon modes, the type and strength of cou- 
pling between the nano-particle and the substrate, and 
finally the initial temperatures of the nano-particle and 
the substrate. A reliable way of studying this problem is 
to carry out state-of-the-art molecular dynamics simula- 
tions by which one can obtain sample specific results. To 
explore the general features of the phononic dissipation, 
however, we propose an Hamiltonian treatment of the 
problem. Since the number of physical ingredients deter- 
mining the dynamics is considerably large, our strategy 



will be to focus on them separately. 

In this work, we present our conclusions concerning the 
dependence of phononic dissipation on internal degrees of 
freedom of the nano-partcile and the substrate, namely 
on the properties of discrete and continuum densities of 
phonon modes of the nano-particle and the substrate. 
Furthermore, we will consider two types of coupling be- 
tween the finite and extended systems. The effect of ini- 
tial temperatures of the parties, besides from the effect of 
temperature difference, is another major topic in its own 
and we leave that discussion to another paper. In the 
present paper, we consider the initial temperatures to be 
zero and limit our attention to the near-equilibrium case 
in the weak coupling regime. Strong coupling regime and 
non-equilibrium cases will also be treated separately. 

The organization of the paper is as follows. In Sec- 
tion HH we describe the physical model. The theoretical 
methods to be used are presented in Section HTT1 These 
are the equation of motion (EoM) technique which in- 
volves Laplace transforms for the solution of the coupled 
differential equations for phonon operators, the Fano- 
Anderson (FA) method which is useful for diagonaliz- 
ing quadratic Hamiltonians, and Green's Function (GF) 
method by which we can incorporate the effect of multi- 
modes into the study. The applications of the given 
methods to different types of coupling and substrates 
having different densities of states will be presented and 
discussed in Section HVl These are mainly the increase 
of decay rate with square of the coupling constant 27 , 
the interaction-specific dependence of decay rate on the 
nano-particle mode frequencies and the effect of neigh- 
boring modes on decay rate of each other. Results ob- 
tained by the above methods will also be discussed in 
connection to classical molecular dynamics (MD) simu- 
lations in Section IIV Al We summarize our conclusions 
in Section [V] 
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Tj. As long as the coupling between the nano-particle 
and the substrate is weak, W^j will be a peaked function 
of c^k and a separate peak will be present around each ujj . 
Depending on the strength of the interaction, the sharp- 
ness of the peaks and the overlap between the neigboring 
peaks will differ. If the coupling is weak enough we may 
neglect the overlaps, namely 

W&Wv - (5) 



FIG. 1: (Color online) A nano-particle with discrete density 
of phonon modes is coupled to a substrate having continuous 
density of modes. 

II. THE MODEL 

We first consider a nano-particle representing a lubri- 
cant molecule or an asperity, which is weakly coupled to 
a substrate (Fig. [T]) . The vibrational modes of the nano- 
particle are excited initially and the excess phonons dis- 
charge to the bulk. The total Hamiltonian of the system 
can be written as 



We assume that the coupling terms of different modes of 
the nano-particle do not overlap, and hence we can treat 
each nano-particle mode separately. 

For the second type of coupling, we consider the cou- 
pling coefficients which scale inversely as the square root 
of the product of the frequencies of the coupled modes, 
i.e. 

Wkj = a {ukUjjY 1/2 . (6) 

The coefficient a stands for the strength of the coupling 
and will depend on the interaction between the nano- 
particle and the substrate. 



H — Hm + Hs + Hms — H + Hms (1) 

where Hm , Hs are the free phononic Hamiltonians of the 
nano-particle and the substrate (or bath) respectively. 
We assume that the harmonic approximation is good 
enough for Hm and Hs, and their spectra are known, 
i.e. 

Hm = hujj afatj, (2) 

3 

H s = 5>co k 6+6 k . (3) 

k 

Here, u)j are the frequencies of the nano-particle modes 
with a,j, being the corresponding annihilation, cre- 
ation operators; Wk are frequencies of the bath modes of 
wavevectors k, and &k , are the corresponding phonon 
annihilation, creation operators. We have omitted the 
constant terms as they do not contribute to the dynamics 
of the system. We consider a single phonon branch with- 
out loosing generality, the formalism can be extended to 
include multiple branches. The interaction Hamiltonian 
Hms is also assumed to be quadratic in phonon opera- 
tors, 

Hms = J2 h ( w V b i a 3 + h - c ) ( 4 ) 

with being the coupling coefficient which is a func- 
tion of Wk and u>j , and has the dimesion of angular fre- 
quency. 

Here we consider two types of coupling. The first one is 
Lorentzian coupling in which the coupling coefficient W^j 
is a Lorentzian with its peak located at u)j and has width 



III. THEORETICAL METHODS 

A. Equation of Motion (EoM) Technique 

The time dependent occupancies of the nano-particle 
modes can be obtained using Heisenberg's equation of 
motion, namely A(t) = i[H,A(t)]/h. The equations of 
motion for the phonon annihilation operators are 

ai{t) = -iwja,(t) - i$^Wki&k(f)i ( 7 ) 

k 

b k (t) = -iwk6 k (i) - iJ2 w kjaj{t), (8) 

3 

that is we have coupled differential equations for each 
operator. Performing Laplace Transformation to both 
equations, a pair of coupled algebraic equations is ob- 
tained 

aiis^s + m) = oi (0) - iJ2 WiMs), (9) 

k 

b k (s) (s + iLo k ) = 6 k (0) - i W kj aj{s), (10) 

3 

where s is the Laplace frequency. Solving for o/(s), one 
obtains 

s + iuji ^ [s + i(dk){s + iwi) 
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Considering the couplings to be non-overlapping (cf. 
Eq. ([5])) we are left with the relation 
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s + iujt + 



s + IUJI+Y, 



(12) 



Having obtained Sj(s) in terms of a/(0) and 6k (0), the 
inverse Laplace transform will yield <i;(t), thus we can 
obtain the time dependent occupancy of the I th mode. 

We convert the summations into integrals over the sub- 
strate modes, and denote them as 
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(13) 
(14) 
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where g{uk) is the phonon density of states for the sub- 
strate, Ii and Ji depend on s, and Ji is an operator. We 
can write Eq. (fT2")) as 



ai(s) 



ai(0) 
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+ itoi + Ii (s) s + iuJi + h (s) 
The inverse transform of a/(s) is 



(15) 



ai(t) 



oi (0) 



2ni J s + iu>i + Ii(s) 

B 



*ds 1 f e st Ji{s)ds 

2w J s + iun + Ii (s) 

B 

(16) 

where the integrals are to be evaluated along the 
Bromwich contour. 

The first and second terms in Eq. (JTBJ) stand for con- 
tributions from the initial excitation of the nano-particlc 
and the initial temperature of the substrate respectively. 
The second term does not contribute to the time depen- 
dent occupations of the nano-particle mode, since the 
initial temperature of the substrate is considered to be 



1. One- and Two-Dimensional Debye Substrates 

For a specific case, we consider the nano-particle to 
have a single mode coupled to a one- or two-dimensional 
Debye substrate which is initially at zero temperature. 
Assuming Lorentzian coupling, namely 



wis - 



we have for Ij{s) 



a 2 T 



1 
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(17) 



FIG. 2: Decay of a single nano-particle mode j coupled to a 
2D-Debye substrate. The coupling is Lorentzian. Occupation 
(nj(t)) at time t is given relative to the initial occupation 
<n;(0)>. 



where d is the dimension of the substrate, cd is the cor- 
responding Debye constant for DOS, gd{u) = Cuu . 
LOj, r are real and positive, and by definition of Laplace 
transformation Re(s) > 0. In the weak coupling regime 
the width of the Lorentzian will be much smaller than 
u>j and u>Dr- so we can approximate the above inte- 
gral by extending the limits of integration to (— oo, oo) 
in which case the integral can be evaluated analyti- 
cally on the complex Wk-plane with the result Ij(s) — 
a 2 c D (ujj - iT/2) d ~ 1 j (s + iu. } + F/2). Performing the 
inverse transformation, one finds 



( nj (t)) = ( nj (0)) 



3 -rt/2 
4IAP" 



(F + A)e At/2 - (r- A)e" 



At/2 



(19) 



where A 2 = F 2 /4 - 4a 2 c D (^- - iT/2) d - 1 . 

Domain of applicability of the EoM technique is quite 
limited due to the fact that inverse Laplace transforma- 
tion is not always possible neither analytically nor nu- 
merically, its advantage is that in certain cases it enables 
us to get analytical results within some approximations. 
Figure [2] shows the decay of a single nano-particle mode 
to a 2D-Debye substrate. 



B. Fano- Anderson (FA) Method 



a 2 Vc D 
2ni 



d-l 



(w k - is)(w k - Uj - i-o)(wk 



Since the Hamiltonian is quadratic in operators, its so- 
7p-(,18) lution is equivalent to diagonalizing a matrix. Exact di- 
''2 ) agonalization of such quadratic Hamiltonians was shown 
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to be possible by Fano2£ and Anderson^ independently 
and the procedure is widely used in atomic physics, solid 
state physics, quantum optics, etc. Here, we will apply 
their method to the problem of phononic dissipation. 

Diagonalization of the Hamiltonian is equivalent to 
finding the dressed operators a(w q ) such that 



H = Hujg a + (ujg ) a {ujg ) . 



(20) 



Substituting Eq. (|21j) into Eq. (|22[) . one ends up with 
a pair of equations 



H(Wq,U}j)(u}q-Uj) = }^ vj^q, Uk) Wkj (25) 
k 

K w q> w k)K-w k ) = y^fijcj^u^W^ (26) 



Since the bare phonon operators a,, 6k form a complete 
set of operators for the combined system, the dressed 
operators a(u>q) can be expanded in terms of the bare 
operators as 

tt(Wq) ^^^Wq.W^Oj + ^ K W q> W k)&k (21) 
3 k 

and they satisfy the eigenoperator equation, 

[a(uJq),H] =?UJqa(Uq). (22) 

Conversely, we find the bare operators by the following 
expressions in terms of the dressed operators as 



a 3 = ^2^*i^i,^j)a(ujq), 
q 

6 k = ^ V*(uJg,UJ^)a{ujq). 



(23) 
(24) 



which can be solved self-consistently to obtain w q . Using 
Eq. (|26p . v can be expressed in terms of jj, as 



v{u)q,UJk) 



UJq - CJ k 



+ S(u>q - Wk)z(Uq) 



^lx{u)q,^)W^ 



(27) 



where P stands for the principal part, and the 5-function 
term accounts for the contribution from the singularity. 
Inserting Eq. (|2"T|) into Eq. (|25p . the following relation 
for fi(ujq,LUj) and z(u>q) is obtained. 



V fi(Uq,U}l)WkjW^ + ^ S(ujq - UJ k )z(u>q)lJ,(ujq,OJl)WkjW^ l 



kZ 



kl 



r 



(28) 



If we consider the nano-particle to have a single mode, 
the relation between /i and z can be written in a much 
simpler form and the dissipation of each mode can be 
treated separately. From this point on, we will use a 
subscript j where necessary denoting that we are working 
on the dynamics of the j mode of the nano-particle. 

Relying on the above reasoning, Zj(ujq) can be ex- 
pressed as 



Zj(u)q) 



aj((Jq) 



<jj(ujq) being the shift in the j th nano-particle mode 

d^k g(u k )\W kj \ 2 



(29) 



r 



LU„ - UJ k 



(30) 



In order to obtain the expansion coefficients [A, v 
phononic commutation relation for a(u)q) is employed. 

[«K),« + K)1 = = S{uJ y^' ] (31) 

gy^q) 



Using the expansion in terms of bare operators (|2ip , and 
Poincare's theorem, i.e. 



P 



Wq — UJk Uqi — U>k UJq' — UJq \ w q~ w k Wq' — ^k 

+ 7T 2 S{uJq - Wk)d(WqJ - UJk), (32) 

the modulus square of n(u>q,ujj) is found as 



\fJ,((jJq,Wj)\ 



K - Uj - (TjK)) 2 + ir 2 g 2 (ujq)\Wq^ 

(33) 

Since the Hamiltonian is diagonal with annihilation, 
creation operators a(u)q), a + (cj q ) and eigenfrequencies 
w q , the time dependence of the dressed annihilation op- 
erator is 

a(w q , t) = Li{uq 7 u j )a j e- iu ' lt + ^ K w q> ^k)&ke" Mq * 
Correspondingly, the time dependence of the nano- 
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particle annihilation operator reads (cf. Eq. ([23]) ) 



Hence the occupancy of the j mode is 



a j (t)= I dwq g(wq) //*(w q , ujj) a(uq) e (35) 



( nj (t)) = ( nj (0)) 



+ / dcj k g(o) k )(n k (0)) 
I 



(36) 



Due to the finite range of substrate DOS g(u>\c), the 
integrals involved in the FA Method are bounded, the 
method allows us to perform calculations for any g(w k ) 
and any type of coupling for a single nano-particle mode. 
In the FA method, the time dependent occupation is 
again separable as contributions from the initial temper- 
ature of the nano-particle and that of the substrate. 



C. Green's Function (GF) Method 



substrate temperature is zero and the phonon modes of 
the nano-particle are empty except for the excitations 
which do not necessarily obey Bose- Einstein distribution. 
That is the initial occupation of a nano-particle mode is 
not a function of temperature. Therefore we make use of 
zero temperature GFs instead of Matsubara formalism, 



d(j,t-t r ) = -i{T t a 3 {t)a+(t')) 



D(k,t-t') 



-i{TA(t)b+(t')) 



(37) 
(38) 



The FA Method is applicable for any coupling type 
and any density of states for the substrate as long as we 
consider a single nano-particle mode. For the multi-mode 
case, the effect of neighboring modes cannot be resolved 
within the above procedures, so we need to develop a 
more generalized method by which we can consider such 
effects. For this purpose we employ GFs. Initially the 



where T t is the time-ordering operator, and the operators 
in Heisenberg picture are distinguished by a hat. 

Since each term in the interaction Hamiltonian in- 
cludes odd number of nano-particle operators, only the 
even order terms contribute in the expansion. The first 
contribution due to the interaction is the second order 
term, 



oo oc 



(39) 



-oo — OO 



whose Fourier Transform gives 

#K) 2 ^^(»'W = d<°>( Wi ) 9 E«( W ,) (40) 

k 

with S (2) (wj) being the second order self-energy. 

1. Nano-F 'article with Single Mode 

First, we need to relate the GF and FA methods, there- 
fore we consider a single nano-particle mode, in which 
the higher order terms can be expressed in terms of the 
second order self-energy and the free GF as 

= d<°>(wi) (l + d (0) S (2) + (#E' 2 ') 2 + ...) . (41) 



For weak coupling, the above series can be written as 

d (0) (^q) 



hence, the retarded function becomes 

1 

d R (Uj,UJq) 



w q - Wj - £( 2 )(o;j-,cj q ) ' 



(42) 



(43) 



The real and imaginary parts of the second order self- 
energy can be separated 



£( 2 >(^-,c q ) = P J 



<2w k #(w k )W? k 

Wq - W k 



iwgiwJW? (44) 



where g(u>) stands for substrate density of states and P is 
for principal part of the integral, and the spectral func- 
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FIG. 3: Diagrams of order In. Solid lines are the phonon 
lines of the nano-particle where the dashed lines are that of 
the substrate, (a) Diagram for the case of single nano-particle 
mode, j. ki stand for the substrate modes, (b) Diagram of 
order 2n when there exists multiple modes (ji) for the nano- 
particle. 



tion is obtained 



A(j, w q ) 



-2ImS( 2 )(j,c Jq ) 



ReE(2)(j, Wq )) 2 + (lmEW(j> q )) 2 ' 



(45) 

We should note here the correspondence between the FA 
and the GF methods. The real part of the second order 
self-energy £( 2 ) is equal to the shift in the j th mode of 
the nano-particle, tjj, obtained previously using the FA 
method, and the square of the imaginary part of E^ 2 ) is 
(ivg(u)q)Wqj) 2 . That is, the FA expansion coefficient fi 
finds its expression in terms of the spectral function as 



|/i(w q ,^)| 2 



2irg(ujq)' 



(46) 



The time dependent GF can be written in terms of the 
spectral function and the time dependent occupancy of 
the j th mode is obtained as 



{nj(t)) = K(0)) 



(47) 



2. Nano-Particle Having Multi-Modes 

As long as the Hamiltonian is quadratic, the primitive 
vertex, from which all diagrams are to be constructed, 
will consist of two phonon lines. That is, each interaction 
point is the intersection of two phonon lines. Since the 
interaction Hamiltonian relates a nano-particle mode to 
a substrate mode only, each vertex contains one nano- 
particle phonon line and a substrate phonon line. So 
the diagram of any order can be constructed (see Figure 
[3Jb). Having obtained the diagrammatic expansion for 
any order 2n, under certain conditions about the coupling 
type, the (2n) th order self-energy term can be expressed 
in terms of the second order term and the free GF of 
the nano-particle modes. If the fraction Wk^/Wk^ is 



independent of ki , the exact form of the self-energy for 
the multi-mode case reads (see Appendix |A"|) . 



w 



(48) 



Once the self-energy is found, the spectral function, 
therefore the time dependent occupancy of the nano- 
particle modes can be calculated. 



IV. DISCUSSIONS AND NUMERICAL 
RESULTS 

In this section, we will apply the above methods to 
Lorentzian and square-root coupling of nano-particles to 
substrates having ID and 2D Debye DOS and discuss 
the dependence of decay rate on properties of coupling, 
and the internal parameters of the nano-particle and the 
substrate. The lowest vibrational frequency of a nano- 
particle increases as its size decreases. In the weak cou- 
pling regime, the width of the spectral function A(j, w q ) 
will be small compared to w q provided that u)j is not close 
to zero, which is already satisfied for nano-particles. In 
this limit, the imaginary part of the self-energy can be 
interpreted as twice the decay rate, -fj = ImS(j, w q )/2. 
Therefore, the dependence of decay rate on the interac- 
tion type and strength as well as on the frequency of the 
nano-particle can be obtained from the spectral function. 

For the case of Lorentzian coupling, the interaction 
strength a is linear with Wkj which shows that the decay 
rate increases with a 2 . If the coupling is a function of 
the separation between interacting atoms of the substrate 
and the nano-particle only, the coupling has the form of 
Eq. ([(j]) with a being proportional to a spring constant 
ki n i connecting the interacting atoms. Since the spectral 
function scales with a 2 , decay rate increases with k 2 nt 
for inverse-square-root coupling case. The k 2 nt law was 
previously obtained using elastic continuum model for 
phononic dissipation at physisorption systems 2 ^. 

The dependence on the nano-particle mode frequency 
is a key issue we wish to emphasize in phononic energy 
dissipation. In Lorentzian coupling case, the decay rate 
is determined by the width of the Lorentzian rather than 
the frequency. On the other hand, for inverse-square-root 
coupling (cf. Eq. ©), it is inversely proportional to the 
nano-particle mode frequency ujj. 

It is evident from equations ([55)) and (|4"5"|) that phonons 
in mode uij decay faster as the substrate DOS at the cen- 
ter of the peak, Wj— Re w q ), increases. A crucial con- 
sequence of dependence on substrate DOS is that, if the 
DOS at the peak of the spectral function tends to zero, 
the spectral function (and |/i(w q , ^j)\ 2 ) has the form of 
a (5-function. In the language of dressed modes, this cor- 
responds to a localized mode, i.e. it does not decay at 
all. Such localized states are also known to occur in e.g. 
solid state physics^ and atomic physics^. For weak cou- 
pling, the real part of the self-energy is small, so the peak 




FIG. 4: (Color online) Effect of neighboring modes. Figures (a), (b), (c) are for lD-Debye DOS, and (d), (e), (f) are for 
2D-Debye DOS with nano-particle vibration frequencies u)i = 0.7u) ma .x, &2 = 0.65 Umax, W3 = 0.55io max , U4, — .45 Umax- 
Figures (a) and (d), (b) and (e) and (c) and (f) show dissipation of phonon occupation for the pairs (lui,lu4,), (0)1,(1)3), and 
(uo\,uj2) respectively. 



of the spectral function is not altered significantly from 
its original position ujj. In other words, lying outside 
the continuum of substrate modes, it is unlikely to be 
shifted into the range where it can decay and vice versa. 
Another important effect about DOS dependence takes 
place when the spectral peak coincides with a van-Hove 
singularity of the substrate DOS, by which the decay rate 
is enhanced abruptly. 

We investigate the effect of a neighboring mode within 
square-root coupling in ID Debye and 2D Debye sub- 
strate densities of states using GF method. We consider 
four nano-particle modes, uj\ = 0.7 u>d, ^2 = 0.65 wd, 
u>3 = 0.55 u>d and W4 = 0.45 u>d- The effect is ana- 
lyzed pairwise, namely, we consider (0)1,0)2), (wi, W3) and 
(u>i, UJ4) as the nano-particle modes keeping other param- 
eters unchanged. That is, we keep o>i constant while 
changing the second mode and investigate dependence 
of decay of tJi-mode on the separation from the second 
nano-particle mode. For both lD-Debye (Figure |4] (a) 
(b) and (c)) and 2D-Debye (Figure H (d), (e), and (f)) 
cases we observe that the decay of excited modes gain a 



retardation as the mode frequencies get closer. A second 
behavior is the enhancement of fluctuations during de- 
cay as the mode frequencies get closer. Both behaviors 
can be understood in terms of the spectral functions. In 
Figure ([Ha) spectral functions of w\— and W3— modes are 
plotted for single-mode (dashed curves) and multi-mode 
(solid curves) GF calculations. It is seen that the overlap 
is negligible and the spectra are not changed consider- 
ably. When the modes are closer (Figure 0b) the single- 
mode spectra (dashed curves) have finite overlap, corre- 
spondingly the multi-mode spectral functions effect each 
other. The Lorentzian shape is distorted and the peak 
of 0J2— mode is enhanced. These result in retardation 
and fluctuations during decay. More precisely, the finite 
overlap of spectra lets the nano-particle to gain phonons 
back which are previously discharged to the substrate. 
This phonon exchange process continues during the dis- 
sipation and gives rise to retardations and fluctuations 
observed in Figure ([4]). 



8 




<°1 






(b) 


w 1 - 

(0 2 

1 
1 

1 

1 


Ml 
III 

M' 
/' I 1 
I J \ 


1 


if ' V 

A* v 

Ik V 

u \ 



0.4 0.6 0.8 

Frequency (in units of w max ) 

FIG. 5: (Color online) Effect of a neighboring mode on the 
spectral function. Dashed curves are the single mode spectral 
functions whereas the solid curves are spectral functions in the 
existence of a neighboring mode, (a) Spectra of u>i = 0.7 ' ui 'max 
and u>3. = 0.55 Wmax for both cases are almost the same, (b) 
Spectra of u)\ = 0.7 uo max and U12 = 0.65 u) max get narrowed 
and distorted when single mode condition is relaxed. 



A. Molecular Dynamics Simulations 

In order to provide another basis on which the results 
obtained from the quantum treatment of the problem can 
be discussed, we refer to a classical treatment of the prob- 
lem. We use a simple but effective approach, where we 
consider the nano-particle/substrate as a cluster/lattice 
of masses and harmonic springs in the first nearest neigh- 
bor approximation, and with the lattice having different 
dimensionalities. The interaction is described by a har- 
monic spring between an atom of the nano-particle and a 
substrate atom. Using the dynamical matrix, the eigen- 
modes of the isolated nano-particle are determined and 
initial energy is loaded to desired modes by giving the 
initial velocities to the atoms in correspondence with the 
modes. In the presence of the interaction between the 
nano-particle and the substrate, the differential equations 
and hence the motion of all atoms are determined in dis- 
crete time steps which are on the order of a femtosecond. 

Since the above stated classical version of the problem 
is not an exact analog to the quantum one, and due to 
quantum vs. classical natures of the two, we compare 
and contrast the basic features of the results emerging 
from them. To be specific, we expect to verify qualita- 
tive issues like dependence on interaction strength, nano- 
particle mode frequency and the effect of neighboring 
modes. 

In agreement with the earlier prediction based on the 
elastic continuum model3£, and the result previously ob- 
tained using the spectral function, the dependence of de- 
cay rate on the interaction strength obeys kf nt law for 
weak coupling. Likewise, the dependence on internal 
degrees of freedom of the nano-particle verifies the pre- 
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FIG. 6: (Color online) MD simulation: ID-substrate and dif- 
ferent nano-particle modes. Higher frequency modes decay 
slower. 



viously obtained result. Keeping the coupling strength 
constant, higher frequency modes decay slower (Figure[(5]). 
We should note that the substrate DOS in the neighbor- 
hood of the nano-particle mode frequency also affects the 
decay rates. Using a ID substrate and choosing the nano- 
particle modes away from the maximum frequency of the 
substrate, the effect of substrate DOS is minimized. Al- 
though the density of substrate phonon modes is higher 
for higher frequencies, the decay rate decreases due to 
nano-particle mode frequency dependence. 

Another property of the dissipation process becomes 
apparent when the dynamics is analyzed when the nano- 
particle has one and two vibrational modes. We con- 
sider a diatomic molecule and a linear triatomic molecule, 
which have one and two vibrational modes along the 
molecular axis, where the interaction is also along the 
molecular axis. The effect of neighboring mode can be 
be analyzed by setting one of the modes of the triatomic 
molecule at the same frequency with the frequency of the 
diatomic one. Exciting only the common frequency of 
both diatomic and triatomic molecules, we compare the 
decay rates keeping the interaction and substrate param- 
eters fixed. In the weak coupling regime, it is observed 
that the decay rate of the common mode does not change 
appreciably. Moreover, exciting both vibrational modes 
of the triatomic molecule does not effect the decay rate 
of the common mode to a great extend. This property 
becomes more apparent as the coupling strength is weak- 
ened. Since the vibrational modes of a triatomic molecule 
along the molecular axis are well-separated this result 
is expected in the light of GF solution of the quantum 
Hamiltonian. The mode localization effect is also tested 
using classical MD simulations. A molecular mode whose 
frequency lies above the maximum frequency of the sub- 
strate has a small but yet finite decay rate. 
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CONCLUSION 
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A theoretical understanding of the phononic dissipa- 
tion from a nano-structure weakly coupled to a substrate 
has been developed using three different methods. The 
EoM technique is able to yield analytical results, but 
has a limited range of applicability because of the fact 
that inverse Laplace transformation is not always possi- 
ble. On the other hand, FA diagonalization is possible 
for any type of substrate density of states and any type 
of coupling, but is restricted to single nano-particle mode 
considerations only. Using GFs, the effect of neighbor- 
ing nano-particle modes can also be investigated. It is 
shown that stronger the coupling, faster is the rate of 
dissipation. Since the width of the spectrum of a single 
nano-particle mode scales with the value of the substrate 
DOS at the shifted frequency of the nano-particle mode, 
we observe that a single nano-particle mode coupled to 
a 2D-Debye substrate decays faster than the one coupled 
to a lD-Debye substrate. This situation can be reversed 
for those frequencies for which ID-DOS is higher than 
the 2D-DOS, namely for low frequencies (larger nano- 
particles). That is, at frequencies at which ID-DOS has 
higher values than 2D-DOS, decay rate of a mode coupled 
to the ID substrate will be higher than that of the mode 
coupled to 2D substrate provided that the remaining fac- 
tors are kept identical. Presence of neighboring nano- 
particle modes effect each other's decay rate when their 
spectral functions have an appreciable overlap. Transi- 
tions between nano-particle modes take place via the sub- 
strate modes, therefore retardation as well as fluctuations 
become important when the modes are close enough. 



This work was supported by The Scientific and Tech- 
nological Research Council of Turkey through Grant No. 
TBAG-104T537. R. T. S. acknowledges financial support 
from TUBA-GEBIP. 



APPENDIX A: SELF-ENERGY FOR 
MULTI-MODES 



When Wkj! / Wk,,7 2 is independent of k, the fourth and 
sixth order contributions can be written as 



£ (4) C/> ? ) 



fes^uuJWSA u; q ))(Al) 



£ (6) (j,- 9 ) = (£rf (0) (ii^)<n J (s (2) 0',- q ))A2) 

By mathematical induction, the (2n) th term is found as 



qj \ h 



x(s (2) 0^ q ))' 



(A3) 
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